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In this addendum to Ref. [I], we give results for the electroweak Sudakov corrections in gauge- 
boson production at large transverse momentum px at proton colliders. In order for the results to be 
easily usable, we provide a simple and accurate parameterization of the corrections as a function of pr 
and the center-of-mass energy y/s. Additionally, we also discuss the dependence of the electroweak 
corrections on the rapidity of the produced boson, and comment on the complications that arise in 
the photon-production case due to isolation requirements. 


In Ref. pT we presented a study of electroweak Sudakov corrections to the cross sections for single W, Z, and 7 
production at large transverse momentum pp, within the framework of Soft Collinear Effective Theory (SCET). Our 
results are based on a factorization formula for the hadronic cross section near the partonic threshold and include 
both Quantum Chromodynamics (QCD) and electroweak corrections. To quantify the effect of including electroweak 
corrections in the cross section, it is convenient to look at the relative difference A a ew of the cross section computed 
with and without electroweak effects, i.e. 
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where a 1 represents the cross section computed including only QCD effects, at a certain order denoted by i , and a l ew 
represents the same cross section including electroweak effects. As shown in Ref. [T], the relative importance of the 
electroweak corrections, given by Aa ew , is essentially independent of the order i which we use to compute the QCD 
effects. This implies that, to good accuracy, one can include electroweak effects in existing pure-QCD computations 
simply as an overall multiplicative prefactor 1 + A a ew . 

The inclusion of electroweak Sudakov effects is essential in order to obtain precise predictions for the pp spectrum 
in the region pp Mw,z, since their effect can be as large as 20% for pp ~ 1 TeV. Therefore, any LHC analysis of 
these processes, as well as studies of future colliders at higher center-of-mass energies, must include them. In Ref. [T) 
we presented the correction factor A a ew for a center-of-mass energy y/s = 7 TeV. Results at yfs = 33 TeV and 
y/s = 100 TeV were presented in Ref. j2j, in the context of the Snowmass study of future colliders. It was found 
that the relative correction A a ew depends only mildly on the center-of-mass energy. For further Run I and future 
Run II LHC analyses it would be useful, to have an expression that allows one to obtain Aa ew for different values 
of y/s. Providing such an expression is the main motivation for this addendum. In order to obtain it, we compared 
the correction factors A u ew at different yfs values and found that their ratios are largely independent of pp. For this 
reason, we first provide a simple parameterization of the results at yfs = 7 TeV, which serves as a reference, and then 
supply the necessary prefactor to convert them to other center-of-mass energies. 

The dependence on pp of A a ew at our reference scale fs/ = 7 TeV can be very accurately parameterized by a 
third-order polynomial, according to 

Aa™(p T ,s 0 ) =a,Q + a\x + afx 2 + afx 3, =: gv(pr), (2) 

where x represents the pp of boson V expressed in TeV, and the af coefficients for the different gauge bosons are 
given in Tab. |TJ both for the central value and the associated error band. The accuracy of the parameterization 
in Eq. ^ is better than 1% for most of the pp values between 100 and 1500 GeV, and the difference between the 
parameterization and the exact result is always much smaller than the theoretical uncertainties of A a ew . Eq. ([2]) can 
therefore be safely used in the whole pp range we study, but should not be employed beyond pp = 1500 GeV. To 
illustrate this, we show in Fig. [T] the exact result for Aa ew (black points) together with the outcome of Eq. |2| (red 
line), as well as the theoretical uncertainty due to scale variations (grey band). 

Next, we provide prefactors fv(fs) that convert Eq. ([2| to other center of mass energies, according to 

A <7™{ P p, s) = gvipr)/fv(Vs) , (3) 


where Aay u (pp, s) represents the difference Aery" at a center-of-mass energy yfs, expressed in TeV. The fv{fs)i for 
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FIG. 1. Aa ew for Z production at yfs = 7 TeV computed using the exact expressions from Ref. [I], black points, and the 
parameterization in Eq. (pi, red line. For reference, we also show the theoretical uncertainty of the result as the grey band (see 
Ref. I] for further details). 


each of the different bosons, are given by: 
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where y = (Vs — yfso)/y/so- Only in the case of W + one observes a slight dependence of the ratio on px and for this 
reason we decided to include also a term linear in px in fyy +; it is understood that px is expressed in TeV. Using the 
fv factors to obtain Aa ew at different energies provides an accuracy better than 1% for most of the px values between 
100 and 1500 GeV, and the difference with the exact result is always much smaller than the theoretical uncertainties. 
We illustrate this in Fig. [2] by comparing ratios of Acr e "’ computed at different energies using the exact result (black 
points) with the fy factors in Eqs. 0-0 (red line). One can also use the fy factors above to convert the upper and 
lower limits of the uncertainty bands at 7 TeV to other yfs values. 

Our finding that A a ew given above has only a mild dependence on fs, is a reflection of the fact that most of the 
electroweak corrections arise in the hard function in the factorization formula, as discussed in Ref. [Tj. The hard 
function evolves from its characteristic scale, which is of order px , to the factorization scale, which is of order Mz, 
and the electroweak corrections generated by this evolution are independent of fs. We note that the predictions of 
our effective-theory setup become less accurate as px —> Mz,w, since the hierarchy of scales on which it is based 
is no longer present. In particular we cannot use it to compute electroweak corrections below px = Mz t w and, 
consequently, one should just set our A a ew factor to 0 for p T < M ZW - Note also that we estimated the uncertainties 
of our results by varying the different scales appearing in the factorization formula independently. When px Mz,w 



upper edge 


V 

a 0 

aX 

aX 

aX 


7.640 

-42.41 

22.88 

-5.609 

7 

4.556 

-22.06 

12.30 

-3.072 

W+ 

9.592 

-46.81 

24.09 

-5.976 

W~ 

9.906 

-49.00 

26.63 

-6.590 



lower edge 


V 

a 0 

a\ 

aX 

aX 


2.302 

-41.49 

22.53 

-5.504 

7 

0.824 

-23.23 

13.44 

-3.428 

W+ 

2.307 

-45.99 

24.97 

-6.244 

w~ 

2.622 

-48.51 

26.84 

-6.637 



central value 


V 

a 0 

aX 

aX 

aX 


2.763 

-40.76 

21.95 

-5.356 

7 

1.713 

-21.68 

12.16 

-3.050 

W+ 

3.816 

-45.25 

23.74 

-5.899 

W~ 

4.074 

-47.60 

25.97 

-6.414 


TABLE I. Values of the af coefficients in Eq. ([ 2 ]), parameterizing the pr dependence of Act™' at fs = 7 TeV, for the different 
gauge bosons V. 
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W* production: ratio 7TeV/33TeV 



FIG. 2. Ratios of A a ew at different energies, for a pp collider, 
over i/s = 33 TeV for Z. The black points are always the exa 
the fv factors in Eqs. 0 - 0 - 
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:t result using the expressions from Ref. pQ, and the red lines are 


s 

L u >> 
b "a 
■o I 

- I % 


0.2 


= 8TeV 


1.0 

100 GeV = p T 

2.2 

100 GeV = p T / 

0.8 

300 GeV = pr 

2.0 

300 GeV = pr - 


'^\\\\\ \ 500GeV = p r 

s 1.8 

500 GeV = p T j 




700 GeV - p T / 1 


\ 700 GeV = pr 

16 

/ / 



V. 

900 GeV = p T / ll 

0.4 

'VY\\\\ \ 900 GeV = p T 

”b 1.4 

1100 GeV = p T / JJ 


\ 1100 GeV = p T 

<1 

1.2 

1300 GeV = pt / ZyS/ 

0.2 

\ 1300 GeV = p T 


1500 GeV = p T 


^g^\\l500 GeV = p T 



on 

.- .. 

n 8 

, 


0.2 


0.4 


y/y m 


0.6 


0.8 


1.0 


y/ym 


FIG. 3. Left panel: Cross section differential in rapidity, normalized to the value at y = 0. Right panel: Size of the electroweak 
corrections as a function of rapidity, normalized to the value of the correction integrated over rapidity. 


the different scales are not really independent and this procedure may result in an overestimate of the errors, but 
since the method itself deteriorates in this limit we recommend to use the error bands given above all the way down 
to pp = Mz. w- 

Our approximate results are for cross sections integrated over the rapdity of the bosons, but to good accuracy the 
corrections can also be used for cross sections with rapidity cuts. To understand this, consider Fig. [3j which shows 
the dependence of the corrections on rapidity for a few fixed pp values. The left panel of the figure shows the cross 
section differential in rapidity, normalized to N = da/dy\ y —o, the value of the cross section at y = 0. The right panel 
shows the size of the electroweak corrections Aa ew (y) as a function of rapidity, normalized to A a ew integrated over 
rapidity. In the horizontal axis, the rapidity is normalized to the maximum value y max , which is given by 


= - In 


S + My - y/ (My - s) 2 - 4 p\s 
Zy/sy/ Pp + M v 


Ii _JL 

2 P 2 t 


( 8 ) 


where My is the mass of the produced boson. Only positive rapidities are shown because the cross section is symmetric 
in rapidity for pp collisions. From the figure we learn that, first of all, the bulk of the cross section comes from the 
central rapidity region, y/y max < 0.8, and, secondly, that the electroweak corrections are essentially flat in this region. 
Therefore, A<j ew is mostly independent of rapidity. Furthermore, the deviations from unity in Figure [3] mostly arise 
for bosons at low pp and therefore correspond to very high rapidities (e.g. at pp = 300 GeV and yfs = 8 TeV, 
2 /max ~ 3.3) which are typically not included in the experimental measurements. The fact that the corrections are 
largely independent of the vector-boson rapidity is demonstrated in Fig. [4j where we have computed them for different 
rapidity ranges relevant in the ATLAS experiment and find that they are almost identical in every case. 
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FIG. 4. Aa ew at 8 TeV for 7 production at different rapidity 
cross section is essentially 0 for pr > 1200 GeV. 
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. In the highest rapidity range shown in the last plot the 


As a final, important point, we briefly want to discuss the effect of experimental cuts beyond the boson rapidity. 
Our computation is based on threshold resummation, i.e. it assumes that the final state only contains the hard 
vector boson, recoiling against a low-mass jet, plus soft gluons and soft photons. In experimental measurements a 
variety of phase-space cuts are imposed and while we hope to capture the bulk of the EW corrections, it is clear that 
the setup of our calculation is not fully realistic. In particular, the EW corrections can change significantly if the 
experimental setup is such that additional, unobserved massive EW bosons can be radiated into the final state since 
their contributions can cancel some of the EW Sudakov logarithms. A second issue related to experimental cuts is 
that near the partonic threshold, the fragmentation contribution in photon production is power suppressed [3], but in 
practice one needs to define what is meant by a final-state photon and depending on the photon isolation requirements, 
a fragmentation contribution will also arise and will involve different EW corrections. The experimental definition 
of the photon final state will also affect which value of electromagnetic coupling should be used when computing the 
the cross section. (The value is obviously irrelevant for the ratio in Eq. 0 we present in this paper, but its choice 
is important when computing the cross section.) In our computation, the appropriate value of the electromagnetic 
coupling is the value at the invariant mass of the vector-boson final state. If the experiment would be able to resolve 
photon jets with arbitrary small mass, the proper value would be p, = 0 and one would use the on-shell, low energy 
value of a e . m . « 1/137, see e.g. the discussion in [3]. However, with the photon isolation requirements at the LHC, 
the photon-jet has a non-zero mass and a value such as a e . m .(Mz) « 1/128 might be more appropriate. 

Acknowledgements: We thank Marat Freytsis, Matt Schwartz and Matt Strassler for discussions. This work is 
supported by the Swiss National Science Foundation (SNF) under the Sinergia grant number CRSII2 141847 1 and 
the US Department of Energy under grant de-sc0013607. The authors thank the high-energy theory group at Harvard 
for hospitality and support. 


[1] T. Becher and X. Garcia i Torino, Phys. Rev. D 88 , no. 1, 013009 (2013) arXiv: 1305.4202 [hep-ph]]. 

[2] K. Mishra et al, arXiv: 1308.1430 [hep-ph]. 

[3] T. Becher and M. D. Schwartz, JHEP 1002, 040 (2010) arXiv:0911.0681 [hep-ph]]. 

[4] A. Czarnecki and W. J. Marciano, Phys. Rev. Lett. 81, 277 (1998) hep-ph/9804252 . 










